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Abstract. 

Context: One of the most debated issues about sub-mjy radio sources, which are responsible for the steepening of the 1.4 
GHz source counts, is the origin of their radio emission. Particularly interesting, from this point of view, is the possibility of 
combining radio spectral index information with other observational properties to assess whether the sources are triggered by 
star formation or nuclear activity. 

Aims: The aim of this work is to study the optical and near infrared properties of a complete sample of 131 radio sources with 
S > 0.4 mjy, observed at both 1.4 and 5 GHz as part of the ATESP radio survey. The availability of multi-wavelength radio 
and optical information is exploited to infer the physical properties of the faint radio population. 

Methods: We use deep multi-colour (UBVRIJK) images, mostly taken in the framework of the ESO Deep Public Survey, to 
optically identify and derive photometric redshifts for the ATESP radio sources. Deep optical coverage and extensive colour 
information are available for 3/4 of the region covered by the radio sample. Typical depths of the images are U ~ 25, B ~ 26, 
V ~ 25 A, R ~ 25.5, / ~ 24.3, 19.5 < Ks < 20.2, J < 22.2. We also add shallower optical imaging and spectroscopy obtained 
previously in order to perform a global analysis of the radio sample. 

Results: Optical/near infrared counterparts are found for ~ 78% (66/85) of the radio sources in the region covered by the deep 
multi-colour imaging, and for 56 of these reliable estimates of the redshift and type are derived. We find that many of the 
sources with flat radio spectra are characterised by high radio-to-optical ratios (R > 1000), typical of classical powerful radio 
galaxies and quasars. Flat-spectrum sources with low R values are preferentially identified with early type galaxies, where the 
radio emission is most probably triggered by low-luminosity active galactic nuclei. Considering both early type galaxies and 
quasars as sources with an active nucleus, such sources largely dominate our sample (78%). Flat-spectrum sources associated 
with early type galaxies are quite compact {d < 10-30 kpc), suggesting core-dominated radio emission. 

Key words. Surveys - Radio continuum: general - Methods: data analysis - Catalogues - Galaxies: general - Galaxies: evolution 



1. Introduction 

The faint (sub-mJy) radio population consists of a mixture of 
different classes of objects. Since the early seventies it has been 
known that the strongest sources are almost exclusively associ- 
ated with either active galactic nuclei (AGNs) or giant elhpti- 
cals, the latter of which are also known as radio galaxies (99% 
above 60 mJy, |Windhorst et al. ||1990| |. More recent work on 
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mJy and sub-mJy sources has revealed that faint sources are 
also found to be associated with normal elliptical, spiral and 
star-forming galaxies, with the early type galaxies being the 
dominant component ( Gruppioni et al. £1999 Georgakakis 



al 



1999[ IMagliocchetti et al. 
Afonso et al. 2006), while at /iJy levels star-forming galaxies 
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prevail (see e.g. Richards et al. [1999 1. 



In spite of the progress made in our understanding of the 
faint radio population, many questions remain open. For exam- 
ple, the relative fractions of the different types of objects are 
still quite uncertain, and our knowledge of their dependence 
on limiting flux density is still incomplete. The reason is, of 
course, that very little is known about the faint ends of the var- 
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ious luminosity functions, and even less is known about the 
cosmological evolution of different kinds of objects. This un- 
certainty is due to the incompleteness of optical identification 
and spectroscopy, since faint radio sources usually have very 
faint optical counterparts. Clearly very deep (/ > 25) optical 
imaging and spectroscopy, for reasonably large deep radio sam- 
ples, are critical if one wants to investigate these radio source 
populations. 

Since the radio emission comes from different types of ob- 
jects an important question is what are the physical processes 
that trigger this emission. It is natural to assume that in the 
case of star-forming galaxies the emission traces the history 
of galaxy formation and subsequent evolution by merging and 
interaction, while the emission in AGNs will reflect black hole 
accretion history. To make matters more complicated, both pro- 
cesses may be present at the same time. 

Although research in this field proceeds slowly due to very 
time-consuming spectroscopy much progress has been made 
in recent years thanks to strong improvement in the photomet- 
ric redshift technique. Several multi-colour/multi-object spec- 
troscopy surveys overlapping deep radio fields have recently 



been undertaken, including the Phoenix Deep Survey ( Hopkins 
erar|[T998| [Georgakakis et al. ||1999[ [Afonso et al. ||2006| l 



and the Australia Telescope ESO Slice Project (ATESP) survey 
( [Prandoni et al. |[2000a|b| [2001j ). In other cases, deep multi- 
colour/multi-wavelength surveys have been complemented by 



deep radio observations (see e.g. the VLA-VIRMOS, Bondi et 
|al. 12003 1 and the COSMOS, [Schinnerer et al."|2006| . 

Multi-frequency radio observations are also important in 
measuring the radio spectral index, which may help to con- 
strain the origin of the radio emission in the faint radio sources. 
This approach is especially meaningful when high resolution 
radio images are available and radio source structures can be 
infeiTed. However, multi-frequency radio information is avail- 
able for very few, and small, sub-mJy radio samples. 

The largest sample with multi-frequency radio coverage 
available so far is a complete sample of 131 radio sources 
with S > 0.4 mJy, extracted from a square degree region ob- 
served at both 1 .4 and 5 GHz as part of the ATESP radio survey 
( [Prandoni et al. |2000a|b|[2006f 

The 1 .4- 5 GHz radio spectral index analysis of the ATESP 
radio sources was presented in the first paper of this series 



( Prandoni et al. 2006 hereafter Paper I). We found a flattening 



of the radio spectra with decreasing radio flux density. At mJy 
levels most sources have steep spectra (a ~ -0.7, assuming 
S ~ v"), typical of synchrotron radiation, while at sub-mJy flux 
densities a composite population is present, with up to ~ 60% 
of the sources showing flat (a > -0.5) spectra and a significant 
fraction (~ 30%) of inverted-spectrum (a > 0) sources. This 
flattening at sub-mJy fluxes confirms previous results based on 
smafler samples (Donnelly et al. |1987 Gruppioni et al. |1997l 
[Ciliegi et al. ||2003j . Flat spectra in radio sources usually indi- 
cate the presence of a self-absorbed nuclear core, but they can 
also be produced on larger scales by theimal emission from 
stars. 

It is possible to combine the spectral index information with 
other observational properties and infer the nature of the faint 
radio population. This is especially important with respect to 



the class of flat/inverted-spectrum sources as it permits us to 
study the physical processes that trigger the radio emission in 
those sources. This kind of analysis needs information about 
the redshifts and types of the galaxies hosting the radio sources. 
A detailed radio/optical study of the sample above is possible, 
thanks to the extensive optical/infrared coverage mostly ob- 



tained in the ESO Deep Public Survey (DPS, [Mignano et al. 
pOOTjlOlsenetal. |2006| ). 

We give a brief discussion of all the data collected so far 
in Sect. |2] followed by a more detailed analysis of the DPS 
optical data in Sect. [3] where we derive the UB VRI colour cat- 
alogue and photometric redshifts for the DPS galaxies in the 
region covered by the ATESP survey, assessing the reliability 
of the photometric redshifts themselves. In Sects. |4] and |5] re- 
spectively, we use the DPS UBVRIJK optical data to identify 
the ATESP radio sources and to derive photometric redshifts. 

A radio/optical analysis of the optically identified radio 
sources is presented in Sect. |6j while in Sect. |7] we discuss the 
nature of the mJy and sub-mJy population on the basis of all 
the radio and optical data available to the ATESP sample. The 
main results are briefly summarised in Sect. [8] 

Throughout this paper we use the ACDM model, with Hq - 
70, Q,„ = 0.3 and Qa = 0.7. 



2. Radio and optical data 

2.1. The ATESP 1.4 and 5 GHz surveys 

As discussed in Paper I, the Australia Telescope Compact 
AiTay (ATCA) was used to image, at 5 GHz, part of the 26° x 1° 
strip of sky previously covered by the 1 .4 GHz sub-mJy ATESP 
survey ( Prandoni et al. |[2000a|b i. In the 2 x 0.5 sq. deg. area 
observed at both 1.4 and 5 GHz a total of 131 distinct radio 
sources are catalogued above a 6cr-threshold (S > 0.4 - 0.5 
mJy) at either 1.4 or 5 GHz (see Table 4 of Paper I). In partic- 
ular we have 89 sources that appear in both the 1 .4 and 5 GHz 
catalogues ( Prandoni et al. ||2000b and Paper I), while the re- 
maining 42 sources are catalogued only at one radio frequency: 
20 sources at 1 .4 GHz and 22 sources at 5 GHz. For the sake 
of the spectral index analysis (see Paper I), we searched for 3cr 
{S > 0.2 mJy) counterparts for these sources at the other ra- 
dio frequency by directly inspecting the (1.4 or 5 GHz) ATESP 
radio mosaics. As a result > 3cr (1.4 or 5 GHz) flux measure- 
ments were provided for 29 additional sources (12 catalogued 
at 1.4 GHz and 17 at 5 GHz), while for 13 (8 catalogued at 
1.4 GHz and 5 at 5 GHz) sources (1.4 or 5GHz) 3cr upper 
limits were estimated. Among the 131 sources catalogued at 
1.4 and/or 5 GHz there are three multiple sources: one is cata- 
logued as triple at 1.4 GHz and as double at 5 GHz and another 
source is catalogued as double at 1 .4 GHz and as a single non- 
Gaussian (extended) source at 5 GHz. 



2.2. The optical/infrared DPS survey 

The 2 X 0.5 sq. deg. area imaged at both 1.4 and 5 GHz as 
part of the ATESP survey overlaps entirely with one sub-region 
(namely the DEEPl sub-region, see below) of the ESO DPS 
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survey. The DPS is a multi-colour survey consisting of both op- 
tical and near-infrared observations. The DPS was carried out 
in the optical {U, B, V, R, I), using the WFI (Wide Field Imager) 
camera mounted at the 2.2m ESO telescope, and in the NIR 
(/, Ks), using the SOFI camera mounted at the ESO NTT tele- 
scope. For a detailed description of the UBVRIJK filters used 



for the DPS we refer to Mignano et al. (2007 1 and Olsen et al. 
|(|2006}. 

The optical (UBVRI) observations cover three distinct 2 x 
0.5sq. deg. regions of sky (named DEEPl, DEEP2, DEEP3). 
Each of the three regions is covered by four 0.5 deg xO.5 deg 
WFI pointings (a, b, c, d). Typical depths of the optical obser- 
vations are Uab ~ 25.7, Bab ~ 25.5, Vab ~ 25.2, Rab ~ 24.8, 
Iab ~ 24.1 ( [Mignano et al. |2007| l. 

The infrared DPS comprises two strategies: shallow K,- 
band (Ksab ^ 21.3) contiguous coverage of about half the WFI 
fields, complemented by deeper J- and A'5-band (Jab ^ 23.4 
and KsAB ^ 22.7) contiguous coverage (4x4 SOFI pointings) 
of the central part of the WFI fields observed in the shallow 
strategy ( .Olsen et al. ||2006 l. In particular for region DEEPl, 
the one of interest for this work, infrared coverage was pro- 
posed for WFI fields DEEPla and DEEPlb in shallow strategy 
(7x7 SOFI pointings), and for the central part of them in deep 
strategy (4x4 SOFI pointings). 

There are some gaps in the optical/NIR imaging of the three 
fields of the region DEEPl (see Table[T]for the summary of the 
observations). Since NIR coverage of each single WFI field is 
obtained with several contiguous SOFI images, the seeing and 
the limiting magnitude values reported in the table for J and 
Ks bands are an average (± standard deviation) of all the SOFI 
images contributing to the WFI field. 

Figure [T] shows the distribution of the infrared SOFI point- 
ings over the WFI fields DEEPla (top) and DEEPlb (bottom) 
for the two available infrared bands (J and K). /T-band cover- 
age is shown in both strategies: shallow and deep imaging (left 
and middle panels). 7-band coverage (right panels) was ob- 
tained only for the deep strategy. The infrared frames are rep- 
resented by the small numbered squares that overlap the corre- 
sponding optical WFI frames (big squares). 

From Table [T] and Fig. [Tjit is clear that DEEPla imaging 
is complete in the optical U, B, and R pass-bands, while no 
imaging is available in the V-band. The /Tj-band imaging, on 
the other hand, is 70% complete in the shallow strategy and 
75% complete in the deep strategy, while 100% completeness 
is reached by the /-band imaging. 

The optical imaging of DEEPlb is complete, except for 
the t/-band imaging, which is slightly shallower than planned 
{miim =~ 24.6). The deep infrared imaging has a good cover- 
age in both filters (>80%), while the shallow /Tj-band imaging 
covers only about 55% of the area. 

It is interesting to note that, even if not complete, the in- 
frared coverage of DEEPla and b is distributed in such a way 
that many of the ATESP radio sources (filled black points in 
Fig. [T]) in the two fields (27 and 26 radio sources for DEEPla 
and b respectively) have infrared information. In particular, 
75% (40/53) of the sources have shallow ^r5-band coverage, 
while deep 7- and /Tj-band infrared data are available for 



Table 1. DPS optical and infrared data status and main at- 
tributes. The table gives in Col. 1 the WFI field, in Col. 2 the 
pass-band, in Col. 3 the seeing, and in Col. 4 the limiting mag- 
nitude (5cr, 2 "aperture, Vega system). 
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Table 2. Main attributes of additional optical imaging obtained 
for fields DEEPla and DEEPlc. Columns as in table [T] 
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24.83 



100% (20/20) of the radio sources located in the central part 
of the fields. 

DEEPlc was only observed in the V band (only down to 
miim ~ 25) and R band, while no observations are available for 
field DEEP Id. 

Reduced images and single pass-band source catalogues 
extracted from both the optical and infrared DPS are described 

( pOOT) and [Olsen et ar| ( [2006j ), re- 



in detail in Mignano et al. 



spectively, and are publicly available at the Centre de Donnees 
astronomiques de Strasbourg (CDS). 

2.3. Additional optical imaging 

Since the DPS was not completed, we have undertaken new 
WFI optical observations in order to collect the missing data 
necessary to have full colour information for region DEEPl, 
and hence for our ATESP radio sources. In this framework we 
have obtained V-band imaging for DEEPla, and U-, B-, I-band 
imaging for DEEPlc. All these new observations were taken in 
collaboration with the group that developed the Garching-Bonn 



Deep Survey (GaBoDS) data reduction pipeUne ( Schirmer et 
[al. |2003[[Erben et al. 12005) and therefore these new data were 
reduced through that pipeline. The main attributes for this addi- 
tional imaging are shown in Table [2] We refer to Hildebrandt et 
[al. [ ( [2006) 1 for a detailed description of the data (both reduced 
images and single pass-band source catalogues). Our multi- 
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Fig. 1. Optical/NIR data coverage for DEEPla (top) and DEEPlb (bottom). From left to right coverage in the K-band (shallow 
and deep strategies) and J-band (deep strategy only). 



colour analysis of ATESP radio sources can rely on full UB VRI 
information for DEEPla, b and c (plus infrared information for 
most of the sources in DEEPla and b). 



2.4. Other optical information 

It is worth mentioning that other optical imaging and/or spec- 
troscopic data are available. The 26 square degree area cov- 
ered by the ATESP survey was chosen to overlap with the re- 



gion where [Vettolani et al. | \\991) made the ESP {ESO Slice 
Project) redshift survey. They performed a photometric and 
spectroscopic study of all galaxies down to bj ~ 19.4. The 
ESP survey yielded 3342 redshifts ( [Vettolani et~aI~|[T998] l, to 
a typical depth of z = 0.1 and a completeness level of 90%. 

In the same region lies the ESO Imaging Survey (EIS) 
Patch A (~ 3° X 1° square degrees, centred at 22''40™, -40°), 
mainly consisting of images in the 1-band, out of which a 
galaxy catalogue (95% complete to / = 22.5) was extracted 
(Nonino et al. 1999[ l. This catalogue allowed us to identify 
~ 57% of the 386 ATESP sources present in that region and 
optical spectroscopy was obtained for a complete magnitude- 



limited (/ < 19) sub-sample of 70 sources (see Prandoni et 
|al. |200ll ). Some VLT/NTT spectroscopy is also available for 
fainter sources in the same region (~ 40 sources with 19 < / < 



21.5, Prandoni et al., in prep.). However, the 3 square degree 
ATESP-EIS sample only overlaps partially with the DEEPl re- 
gion, covering the fields DEEPlc and DEEPld. 

This paper mainly focuses on the radio/optical analysis of 
the 85 ATESP radio sources located in DEEPla, b and c, for 
which deep multi-colour optical/NlR information can be ex- 
ploited. However, whenever considered useful, we include in 
our discussion any optical data (imaging and/or spectroscopy) 
available to the ATESP sources located in DEEPld. Such data 
may come either from the observations mentioned above, or 
from the literature. 



3. Multi-colour analysis of DEEP1 DPS data 

A general discussion of the DPS optical imaging is provided in 
Mignano et ar| ( |2007| l and in |Hildebrandt et al. | ( p006l l, where 
the global quality of the photometry obtained through the EIS 
and the GaBoDS pipelines, respectively, is discussed. Here, we 
focus our attention on region DEEPl, which is the region of in- 
terest of this work. A careful analysis of the photometry of the 
single pass-band images covering region DEEPl is very im- 
portant since we will use this data later to estimate photometric 
redshifts for the ATESP radio sources. Also very important is 
the recipe followed to produce the optical colour catalogues. 
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since reliable galaxy colours are crucial to get reliable photo- 
metric redshifts. 



3.1. Colour catalogues 

We used the available UBVRI images to derive overall optical 
colour catalogues for DEEP la, b and c. 

To obtain a good quality colour catalogue it is clear that 
one should use UBVRI images reduced in a consistent way. 
The optical images available to this work were reduced with 
different pipelines: the EIS pipeline for the images obtained in 
the framework of the DPS survey and the GaBoDS pipeline for 
the images obtained later on (UBI images for DEEPlc and V 
images for DEEP la). In order to avoid internal inconsistency, 
we therefore decided to refer to the EIS reduction for DEEP lb 
(see 'Mignan o et al. |[2007 1 and to the GaBoDS reduction for 
both DEEPla and DEEPlc (see |Hildebrandt et al. |2006 ). 

The technique of reference imaging (see below) was 
adopted in this work to produce the colour catalogue, since it 
provides the most accurate colour estimates, through the mea- 
suring of the source flux within the same area in any of the 
dififerent pass-band images. This is especially important for ex- 
tended objects. 

We selected as the reference image the best seeing single 
pass-band image. Such a choice allows us to minimise the ef- 
fect of very close pairs of objects, which are not resolved due to 
poor seeing. The /-band and R-band images for field DEEPla 
have very similar seeing values and the colour catalogue of 
DEEPla was extracted by using the R-hwd image as the ref- 
erence. For DEEPlb and DEEPlc best seeing is measured for 
I- and R-band images, respectively, and the choice of the ref- 
erence image was done accordingly. 

We run SExtractor (ver. 2.3, [Bertin & Arnouts |[T996l l in 
the so-called double image mode: detection and object aper- 
tures were based on the reference image, followed by isophotal 
magnitudes measured in the same aperture for each detected 
object on the other pass-band images separately. 

The optical colour catalogues were then cross-correlated 
with single pass-band catalogues extracted from the J and 
images, whenever available. Since the infrared catalogues over- 
lap, it may happen that the same optical object is identified in 
more than one infrared catalogue. In such cases, the infrared 
object with the lowest magnitude eiTor was selected. 

We did not include the NIR information in the colour cat- 
alogue production from the beginning since a) it is available 
only for limited sub-regions of DEEPla and DEEPlb fields 
and b) the data are taken with a different instrument and tele- 
scope (SOFI at the 3.6 m) and reduced through a specific EIS 
pipeline. 



3.2. Photometric redshifts 

The success of photometric redshift estimate routines strongly 
depends on the accuracy of the photometric calibration in the 
various pass-bands and on the accuracy of the colour estima- 
tion. In Mig nano et al. IpOOTj l and [Hildebrandt et al. | ( |2006| l, 
comparisons between the UBVRI colours of stars in the vari- 
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Fig. 2. U-B vs. B-V colour diagram for stars in field DEEPlb. 
No coiTection applied (top), correction Ucorr = -0.15 applied 
(bottom). Green points refer to DPS stars, black dots to mod- 
elled stars. 



ous regions covered by the DPS, and the ones expected from 

2005[ l, were presented to 



a theoretical model (Girardi et al 



check for the presence of possible systematic offsets. Here, we 
report on the results obtained specifically for the DEEPl re- 
gion, which is the one of interest to our radio/optical study. 
From the colour-colour diagram analysis very good agreement 
was found between the catalogued star colours and the theo- 
retical expectations, except in the case of the U-band for field 
DEEPlb, where an off'set of ~ 0.15 mag is present (see Fig.|2j 
top panel). After correcting for this offset, a good overlap be- 
tween observed and expected colours is obtained (see Fig. |2j 
bottom panel). 

We have also checked the optical (WFI)-infrared (SOFI) 
colours of the DPS stars, and no appreciable off'set was seen. 
This is shown in Fig.|3] where V - R vs. R - J and R - are 
plotted for DEEPlb, chosen as reference. 

We also used any spectroscopic data available from the lit- 
erature in this region to analyse the impact of both the correc- 
tion applied in the U-band for DEEPlb and the use of NIR 
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Fig. 3. Optical vs. infrared colour-colour diagram for stars in 
DEEPlb: V-Rvs.R-J (top) and V - 7? vs. R - K, (bottom). 
Green points refer to DPS stars, black dots to modelled stars. 



colours (when available) in the determination of galaxy photo- 
metric redshifts. 

It is important to note that most spectra come from the ESP 



redshift survey (see Sect. 2.4 1, which covers a limited redshift 



range (z < 0.3, Vettolani et al. [1998 1 and therefore the present 
comparison mainly probes the most local galaxies of the DPS. 

Photometric redshifts for 88 galaxies with spectroscopy in- 
formation present in fields DEEP la, b and c were estimated 



using the public photometric redshift code Hyperz (Bolzonella 



et al. 20001, by using both the templates created from the 



synthetic stellar libraries of Bruzual & Chariot ( 1993[l, here- 
after EC, and the empirical ones compiled by Coleman et 



|al. I ( |1980[ ) to represent the local galaxy population (hereafter 
CWW). We stress that such galaxies are not necessarily associ- 
ated to ATESP sources. 

From this analysis we found a clear improvement in the 
photometric redshift determination when correcting for the sys- 
tematic off'set in U-band photometry. A further improvement is 
obtained when adding the NIR {J and K) information (when 
available) to the UBVRI colour catalogue. The Zpiwi vs. Zspec 



linear fit slope gets closer to unity (a - 0.93 + 0.05) and the ob- 
ject distribution around the Zphot - Zspec line gets narrower. The 
final Zphot — Zspec diagram is shown in Fig.|4] Dotted lines indi- 
cate the range that contains 95% of the objects (Zp/,„, - Zspec± 
0.1). Such a range, albeit rather large, is adequate for this kind 
of study, where errors in luminosity determinations of AlogL 
of the order of ^ 0.5 are acceptable. The horizontal error bars 
are not shown since they are negligible: ESP redshifts are char- 
acterised by errors of the order of ~ 60 km/s on the measu red 
recession velocity, i.e. Az ~ 2 ■ 10 ( Vettolani et al. |1998 i. 

As a final remark, we stress that photometric redshifts 
shown in Fig. |4] were obtained using different template sets 
for diff'erent redshift ranges: galaxies with spectroscopic red- 
shift < 0.1 were fitted by CWW templates, while objects at 
Zspec ^ 0.1 by BC templates. This choice, as expected, turned 
out to provide the best redshift estimates over the two redshift 
ranges. 

Figure |5] shows the galaxy photometric redshift distribution 
obtained from the optical UBVRI colour catalogues in field 
DEEPlb, chosen as reference. Most of the galaxies lie at z < 1, 
as expected, with a significant number of objects extending up 
to z ~ 3. On the other hand, the excess at z ~ 5.5 is mainly due 
to objects classified by Hyperz as Sc galaxies and is clearly 
spurious. For such objects the photometric redshift determi- 
nation clearly fails. Noteworthy are the two naiTow peaks at 
z ~ 0.7 and z ~ 1.5. The latter is also present in fields DEEPla 
and DEEPlc, and most probably indicates a degeneracy in the 
Hyperz routine, due to the fact that the spectral range covered 
by UBVRI-bands at z > 1.5 does not probe the Balmer 4000 A 
break. More interesting is the peak at z ~ 0.7, mainly composed 
of early type galaxies, which is not replicated in DEEPla and 
c, possibly indicating the presence of real large scale structure. 



4. Optical identification of the ATESP radio 
sources 

In the following we present the cross-correlation between the 
ATESP radio sources in fields DEEPla, b and c and the multi- 
colour optical/NIR catalogues described in Sect. 



3.1 



In the literature diff'erent statistical techniques are used to 
cross-correlate radio and optical catalogues, from the simplest, 
distance-only based criterion, which considers as good any 
identification within a certain fixed radio-optical distance, to 
more sophisticated techniques, like the likelihood ratio crite- 
ria, based on the probability that a given source, at a certain 
distance and with a certain magnitude, is the true optical coun- 
terpart of the radio source (e.g.'de Ruiter et al. '1977''CiUegi et 
al. |2003t [Sullivan et ar[2004, Simpson et al. .2006.) . 



In Mignano et al. ( 2007 1 a prehminary optical identifica- 
tion of the ATESP sources with the DPS catalogues was pro- 
posed, based on distance alone. However, while this choice 
proves to be appropriate for shallower optical databases (see 
e.g. the ATESP-EIS case, Prandoni et al. \2QQT I, it is not very 



reliable when dealing, like in this case, with very deep (and 
therefore crowded) optical catalogues. Hence, we adopt the 



likelihood ratio technique in the form described by Sutherland 
& Saunders \^T992^ and |Cihegi et"!!"! ( [20031 ). 
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where C is the completeness and R the overall reliability of the 
sample ( de Ruiter et al. 1 1977 1. 



Fig. 4. Zpi,ot vs. Zspec for the 88 galaxies in DEEPla, b, c fields 
with spectra available. f/-band magnitudes of DEEPlb objects 
are corrected and NIR information is used, when available. The 
solid line indicates Zphot - Zspec and the dotted lines indicate the 
range that contains 95% of the objects. The error bars represent 
the limits of the confidence intervals at 68% 




Fig. 5. Galaxy photometric redshift distribution for field 
DEEPlb. 



The likelihood ratio LR is defined as the ratio between the 
probability that the source is the correct identification and the 
coiTesponding probability that the source is a background, un- 
related object. A threshold value L,/, of the Ukelihood ratio is 
assumed, above which a counterpart is considered as a good 
identification and below which is dismissed as spurious. 

The sample of accepted identifications thus consists of all 
the radio-optical associations that have LR > L,/,. L,i, was cho- 
sen to be the value of LR that maximises the function (C+R)/2, 



4.1. Optical identifications 

The ATESP radio sources were identified in the same reference 
pass-band (R or I) as chosen to derive the colour catalogues of 
DEEPla, DEEPlb and DEEPlc. 

Before proceeding with the optical identification, the pres- 
ence of possible systematic offsets between the radio and the 
optical astrometry was verified. We note that radio positions 
always refer to 5 GHz catalogue positions, unless the source is 
catalogued only at 1.4 GHz (i.e. S'peakCS GHz)< 6cr, see Paper 
I), while optical positions refer to the reference pass-band cat- 
alogue. 

( 2QQ7\ , where a preliminary 



As shown in Mignano et al. 



analysis was given, the median radio-optical offsets for our 
sample are < A RA>= -0. "213 and < A Dec>= -0. "073. The 
source radio positions were coiTected for such median off'sets 
before proceeding with the radio source optical identification. 

In computing the LR value for each optical counterpart, the 
radio and optical positional uncertainties have to be taken into 
account. Here, we adopted Icr positional errors appropriate for 
the ATESP (see [Prandoni et al. ||2000b| ) and for the DPS (see 
[Mignano et al. |[2()07 1 catalogues. In addition, we need to as- 
sume an expected a priori identification rate (Q). We adopted 
Q - 0.7, i.e. 70% of the radio sources are assumed to be truly 
identified down to the limiting magnitude of the optical cata- 
logues. This choice is based on previous radio-optical identifi- 
cation studies undertaken down to similar optical depths (see 
e.g. Ciliegi et al. |2005| [Sullivan et al. |2"004| . 

Figure |6] shows the distribution of LR values as a function 
of radio-optical offsets for the 85 radio sources in DEEPla, b, 
and c. As expected, LR decreases going to large radio-optical 
offsets (> 1 "), and the identifications become less reliable. 
The horizontal solid line represents the assumed threshold LR 
value, above which optical counterparts are considered as good 
identifications. The adopted threshold value, L,/, = 0.3, was 
chosen in agreement with similar works reported in literature 
(e.g. Ciliegi et al. |2005) l. It is worth noting, however, that most 
of the sources have LR values » 10 (see Fig.[6]l, which means 
that most of the optical identifications have very high probabil- 
ity of being real. As reported in Table [3] 60 radio sources in 
DEEPla, b and c were identified down to L,i, = 0.3, (see Col. 
3). 

In order to check the robustness of this identification tech- 
nique and its dependence on the assumed parameters, the like- 
lihood ratio analysis was repeated using different values of Q in 
the range 0.5-1.0. No substantial difference in the final number 
of identifications and in the associated reliability was found. 

The contamination due to possible spurious identifications 
with LR > Lfi, was estimated by shifting the coordinates of the 
radio sources by several random offsets and then repeating the 
identification procedure. The average contamination rate i%sp) 
was 7.4%, 6.8% and 6.3% for DEEPla, DEEPlb and DEEPlc, 
respectively (see Table [3] Col. 5). 
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Fig. 7. The three extended radio sources included a posteriori in the identification sample: two with double radio morphology, 
ATESP5 J225034^01936 (left panel) and ATESP5 J225426^02442 (middle panel); and one wide angle tail source, ATESP5 
J225505-401301 (right panel). Grayscale: optical image in the reference band (R for DEEPla and / for DEEPlb). Contours: 5 
GHz flux density. 
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Fig. 8. The three sources identified in bands other than the reference one: ATESP5 J224827-402515 (left panel), identified in 
/-band; ATESP5 J225443-401147 (middle panel), identified in /T.-band; and ATESP5 J22551 1-4015 13 (right panel), identified 
in 7-band. Grey scale: optical or NIR image. Contours: 5 GHz flux density. 



Table 3. The sample identification statistics. Col. 1 lists the 
field name. Col. 2 the number of radio sources. Col. 3 the num- 
ber of objects identified using the LR technique. Col. 4 the com- 
pleteness of the sample. Col. 5 the average contamination rate. 
Col. 6 the number of additional identifications. Col. 7 the total 
number of identifications, and Col. 8 the identification rate. 



Field 






C 




id 


N- 


%,d 


DEEPla 


27 


16 


98.6 


7.4 


4 


20 


74.1 


DEEPlb 


26 


21 


99.1 


6.8 


1 


22 


84.6 


DEEPlc 


32 


23 


99.0 


6.3 


1 


24 


75.0 




85 


60 


98.9 


6.8 


6 


66 


77.6 



4.2. Additional identifications 

The optical counterparts of the ATESP radio sources were 
all visually inspected on the corresponding optical reference- 
band images, giving particular attention to the multiple and 
non-Gaussian radio sources, where radio positions might not 
precisely coincide with the host galaxy core. From this in- 



spection we decided to include among the identified radio 
sources three extended radio sources with LR values lower 
than the threshold: in all cases, the optical counterpart lo- 
cated close to the radio barycenter is very likely to be the host 
galaxy of the radio source. Two (ATESP5 J225034-401936 and 
ATESP5 J225426-402442) are classical double radio sources 
(see Fig.|7| left and middle panels) and one (ATESP5 J225505- 
401301) has the morphology of a low surface brightness wide 
angle tail (WAT) source, i.e. an extended radio source located 
in a cluster (see Fig. [7] right panel). This hypothesis is sup- 
ported by the fact that ATESP5 J225505-401301 is located in a 
crowded optical field with several optical galaxies in the field 
having similar photometric redshifts (see Sect.|3]for the deriva- 
tion of photometric redshifts for the optical sample). 

In addition, we checked for any possible additional iden- 
tification in pass-bands other than the reference one {R or /). 
We found that only one extra identification (source ATESP5 
J224827-402515 in DEEPlc field) could be recovered when 
the reference optical catalogue was extracted from the /-band 



A. Mignano et al: The ATESP 5 GHz radio survey. II. 



9 



• • • . - 


. • % 

• * • 






■ J 

• 




-1 
1 



Table 4. Identification rate in our sample and in other deep 
radio fields: VVDS-VLA ( jCiUegi et al"~|[2005]l, Phoenix sur- 



vey ([Sullivan e t al. |2004|l, VLA-LH ( [Ciliegi et al. |2003l ) and 
ATESP-EIS ( ,Prandoni et al. |[200T] ). Col. 1 gives the sample 
name, Col. 2 the radio flux limit, Col. 3 the number of radio 
sources present in the sample. Col. 4 the area covered by the 
radio-optical data. Col. 5 the limiting magnitude (/), and Col. 6 
the identification rate. 



r (arcsee) 



Survey 






Area 


I/im 


%,d 




(mJy) 




(sq.degr.) 






ATESP-DEEPl 


0.4 


85 


0.75 


24.3 


77.6 


VVDS-VLA 


0.08 


1054 


I 


24.5 


74.0 


Phoenix 


O.I 


839 


3 


24.5 


79.0 


VLA-LH 


0.05 


63 


0.03 


24 


92.0 


ATESP-EIS 


0.4 


386 


3 


22.5 


57.3 



Fig. 6. LR values vs. radio-optical offsets for the 85 radio 
sources located in DEEP la, b and c. The horizontal solid line 
indicate the LR value L,/, - 0.3, above which counterparts are 
considered as good identifications (see text for details). 



image (see Fig.jSJ left panel). This identification was originally 
missed due to the fact that a larger region around the bright 
star close to the object was masked in the /?-band image. This 
means that in general the reference images were chosen appro- 
priately for our scientific application. 

A similar check was performed for the J- and Ks-hand in- 
frared images, available for DEEP la and b. For two sources 
(ATESP5 J225511-401513 and ATESP5 J225443-401147), 
with no optical counterpart in the DEEPla optical images, a 
possible counterpart was found within a distance of 2 " in the 
infrared Ks- or /-band images (see Fig. [8] middle and right 
panels). These objects have extremely red colours {R-Ks > 5), 
probably caused by either high redshifts or reddening due to 
dust. 

Including the six objects discussed above (Nf^'' in Table|3j, 
the final identification sample is composed of 66 objects, cor- 
responding to an identification rate of 77.6% (see last column 
of Table [3]l. On average the completeness C is 98.9% and the 
contamination rate is 6.8%. Both these quantities refer to the 
sample of 60 identifications, statistically defined on the base of 
the likelihood ratio technique. A summary of the sample iden- 
tification statistics is given in Table [3] 

A list of all the identified radio sources is given in Table |6] 
The six objects included a posteriori (see discussion above) are 
added at the bottom of the table. 

A comparison with other similar radio/optical studies is 
shown in Table [4] It is worth noting that the identification rate 
of our sample is consistent with the ones found in similar radio- 
optical studies taken from the literature. Of particular interest 
is the comparison with the identification rates reported for the 
VVDS -VLA sample ([Ciliegi et al. |2005] l, and for the Phoenix 
survey ( |SuUivan et al. [2004 ^, where the radio/optical analysis 
was performed down to the same optical depth. 



It is also interesting to compare the present study with the 
shallower ATESP-EIS sample, where optical identifications 
were searched down to / = 22.5 (see [Prandoni et aT71|2001[ ). 
The identification rate increases from ~ 57% of the ATESP- 
EIS to 78% of the ATESP-DEEPl, demonstrating the need for 
deep follow-up surveys to properly identify the mJy/sub-mJy 
radio population. 

5. Photometric redshifts for the identified ATESP 
radio sources 

The optical spectra of radio sources are not exhaustively repre- 
sented by the standard "stellar" templates used for normal in- 
active galaxies (Ellipticals, SO, Spirals, Irregulars, Star-bursts). 
We therefore added to the standard template spectra provided 
by Hyperz (BC and CWW) a set of spectral templates derived 
from the SDSS (Sloan Digital Sky Survey, |York et al. ||2000| l 
quasar samples. In particular we added templates for: a) blue 
quasars (QSO, Hatziminaoglou et al. 2000); b) composite red 
quasars (REDQ, Richards et al. ,2003) ; and c) composite broad 
absorption lines (BAL) quasars (BALQ, Reichard et al. [2003 1. 
These templates are available at the SDSS web pages. 

For each identified ATESP radio source Hyperz was used to 
provide a possible Zphoi (and corresponding reduced proba- 
bility) for each set of templates (CWW, BC, blue, red and BAL 
quasars). Then, the "best" (highest probability) Zp/,o; was se- 
lected as the correct one, together with the coiTesponding spec- 
tral type. 

In eleven cases it was not possible to assign a reliable 
Zphoi and spectral type to the optical radio source counter- 
part. Typically these are very faint objects (mag>24 in detected 
bands), or objects with very limited colour information (e.g. 
detected only in NIR pass-bands), or objects with bad pho- 
tometry due to nearby star and/or deblending problems. One 
of these cases could be recovered thanks to the availability of 
spectroscopy information (source ATESP J224958-395855). 

In summary, it was possible to assign a redshift and a spec- 
tral type to 56 of the 66 radio sources identified in DEEPla, 
b and c (85%). However, if we restrict our analysis to a 
magnitude-limited / < 23.5 complete sample, we get a success 
rate of 97% (56/58 objects with redshift determination). 
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The relevant spectral parameters obtained for the 56 radio 
sources for which a redshift and type could be assigned, are 
reported in Table |6] 

Spectral types reported in Table |6]( Col. 13) are defined as 
in |Prandoni et al. | ( p50T] l): 

1. Early type spectra (ETS): ellipticals, early spirals (bulge- 
dominated Sa); 

2. Late type spectra (LTS): late spirals (Sb, Sc, Sd) and irreg- 
ular Magellanic (Im) galaxies; 

3. SB: star-burst galaxy spectra (typical of HII regions); 

4. PSB: post star-burst galaxy spectra (K-hA and E-hA galax- 
ies); 

5. AGNs: objects with evident characteristics of either Seyfert 
1, Seyfert 2, or quasar spectra (respectively labeled as Syl, 
Sy2, Q); 

For the 14 objects with optical spectroscopy available, we 
have in general a good match between spectral and photomet- 



Table 5. The ATESP sample composition. 



ric redshifts (Az ^ 0.1, see Sect. 3.2 1. The exception is ATESP 
J224803-400513, for which we have a Az » 0.1. In this 
case, the obtained photometric redshift {Zphor = 1-0) is much 
lower than the spectroscopic value of 1 .72 (the published value, 
Zspec = 2.33 by Prandoni et al. |200T was over-estimated). 

We also find very good agreement between photometric and 
spectral types. There is only one case (source ATESP J225400- 
402204) where the photometric type (very old Sa) disagrees 
with the spectral type (LTS). However, note that passively 
evolving single bursts (defined as Burst in Hyperz BC tem- 
plates) can be considered early or late, depending on their age. 
As a general rule, very old galaxies (age > 1 Gyr) are included 
among the ETS, while very young galaxies (age < 0.1 Gyr) 
are included among the LTS. For intermediate cases (ages be- 
tween 0. 1 and 1 Gyr) the classification is not straightforward, 
from wide-band information on the continuum shape alone. It 
is difficult to distinguish between LTS, ETS and PSB, with no 
information on the presence of narrow absorption and/or emis- 
sion lines. For the sake of simplicity we decided to make a 
sharp separation between LTS and ETS at age = 0.3 Gyr, with 
the caveat that among such objects we could have some mis- 
classification. The value of 0.3 Gyr was chosen from a compari- 
son between spectral type and Burst age in the few cases where 
spectroscopy was available. One probably mis-classified ob- 
ject is source ATESP5 J22532 1-4023 17 (Burst age 0.18 Gyr), 
which has linear size ~ 200 kpc and extended radio morphol- 
ogy, clearly indicating an AGN origin of its radio emission (see 
Fig. [13] middle panel). 

The final classification given to the optical counterpart on 
the base of the present discussion is reported in last column 
of Table [6] The classes are defined following the spectral type 
definitions listed above. In the one case where photometric and 
spectral types disagree, we rely on the latter to define the object 
class. 

To further check the reliability of our photometric redshift 
determinations, we compared our redshift distribution with the 
one expected for ETS on the base of the well-known K-z cor- 
relation found for radio source host galaxies (see e.g. ' Willott| 
et al. |20()3] l and with the R-z relation found for host galaxies 



Sample 






ETS 


LTS+SB 


AGNs 


UNCL 








(%) 


(%) 


(%) 


(%) 


ATESP- 


-EIS 


19 


49 ±8 


43 ±8 


9±3 




ATESP- 


-DEEPl 


23.5 


64± 10 


19 ±6 


14 ±5 


3±2 



of gigahertz peaked spectrum (GPS) radio sources (see Snellen 
et al. 1 1996) |Rigby et al. |2007| ). We find that the photometric 
redshift distribution obtained for objects classified as ETS in 
our sample follows, within a Az ~ 0. 1 - 0.2 dispersion, the one 
expected on the base of the quoted relations. 

6. The ATESP-DEEP1 source properties 

We exploited the photometric redshift and spectral type de- 
terminations for the ATESP sources in the DEEPl a, b and c 
regions, to study the composition of the ATESP sample and 
the radio/optical properties of mJy and sub-mJy sources. It 
is important to note that the ATESP-DEEPl sample overlaps 
with the ATESP-EIS sample. This means that photometric (or 
spectroscopic) redshifts obtained for sources in DEEP la, b and 
c could be in principle complemented by the sparse spectro- 
scopic information available from the ATESP-EIS for DEEP Id 
sources. Nevertheless, in the following it was preferred to limit 
our analysis to the sources in DEEP la, b, and c, (see Table |6]), 
which represent a much more reliable sample, thanks to the 
very high identification/redshift determination statistics. 

Of the 66 identified radio sources, it was found that 37 are 
ETS, 8 are quasars (AGNs), 10 are LTS, 1 is a SB, while 10 
objects could not be classified (UNCL). 

In Table [5] the ATESP-DEEPl composition is compared 
with the one found for the "brighter" ATESP-EIS sample (70 
objects with complete spectroscopy down to / = 19, Prandoni 
pOOT). The ATESP-DEEPl sample provides insight into 



etal. 



the composition of the faint radio population associated with 
optically faint galaxies, albeit in this comparison we restrict 
our analysis to the magnitude-limited sample of 58 objects with 
/ < 23.5, to reduce the number of unclassified obje cts. Table[5 
shows that, as suggested by previous studies (e.g. Gruppioni 
let al. II 1999) |Prandoni et al. ||200l 1), the contribution of star- 
forming (LTS plus SB) galaxies decreases dramatically with 
the magnitude, going from 43% of the ATESP-EIS "bright" 
(/ < 19) sample to 19% of the "deeper" (/ < 23.5) ATESP- 
DEEPl sample. The fraction of ETS and AGNs, on the other 
hand, increases going to deeper magnitudes, even though the 
statistical uncertainties are large. 

6.1. Redshift distribution 

Figure [9] shows the redshift distribution of the 56 ATESP ra- 
dio sources in regions DEEP la, b and c, for which a reliable 
redshift estimate was obtained (see Table |6|l. Whenever spec- 
troscopy is available (14 objects), we rely on the spectral red- 
shift determination. The distribution of ETS shows a significant 
peak at z=0.4, with a tail extending up to z~ 2, while, as ex- 
pected, quasars have typically higher redshifts 1 < z < 2, and 
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Fig. 9. Redshift distribution for the 56 radio sources in the 
ATESP-DEEPl sample with photometric redshift determina- 
tion. The sample is divided into four different classes. From 
top to bottom: AGNs, ETS, LTS, and star-burst galaxies. 



Fig. 10. 1.4 GHz radio power distribution for the 56 ATESP- 
DEEPl radio sources with photometric redshift determination. 
The sample is divided into four different classes. From top to 
bottom: AGNs, ETS, LTS and star-burst galaxies. Light shad- 
ing indicates two upper limits. 



LTS are found at z«l. This reflects the fact that radio sources 
triggered by star formation are usually characterised by lower 
radio powers than sources triggered by AGN activity (see also 
Fig. 10 1. The only star-burst galaxy in the sample has either a 
redshift (z ~ 2) or a radio power (Pj 4 ghz close to 10^^ W/Hz, 
see Fig{TO| that is much higher than for the LTS galaxy pop- 
ulation. While this could be due to evolutionary effects in the 
population of the radio-selected star forming galaxies, it is also 
possible that the photometric classification is wrong. In fact the 
SB spectra are notoriously similar to narrow-line AGN spectra 
(Seyfert 2) and photometric techniques based on wide-band 
colours could easily fail in classifying such objects. In addi- 
tion, the photometric routine applied to this sample (Hyperz) 
does not provide template spectra for Seyfert 2 galaxies. 

6.2. Radio and optical luminosities 

For the 56 ATESP radio sources in DEEPla, b and c with red- 
shift determination we derived radio and optical/NIR luminosi- 
ties. Radio powers were K-corrected by using the 1.4-5 GHz 
radio spectral index of each source (see Table [6]l, while abso- 
lute magnitudes (computed by Hyperz) were K-corrected on 
the base of the optical spectral type. 



Figure 10 shows the 1.4 GHz radio power distribution for 
the sample. Again, the four classes (AGNs, ETS, LTS, and 
SB) are shown separately. ETS galaxies mostly have 23 < 
logf (W/Hz) < 25, which are typical values of FRI radio 



sources (Fanaroff & Riley 1974 1, while AGNs are, as expected, 
characterised by higher radio powers (10^^ - 10^* W/Hz). LTS 
galaxies, on the other hand, have low radio powers, with 7/10 
having P < 10^* W/Hz, typical of radio sources triggered by 





10 - 



I ' ' ' I ' ' ' I 

AGN 



-|— I — I I I I — I — H 




LTS 



I I I I I I I 
SB 



Star formation (see e.g. Condon & Broderick ,1988) 1. 



Fig. 11. Absolute Z?-band magnitude distribution. The sam- 
ple is divided into four different classes. From top to bottom: 
AGNs, ETS, LTS and star-burst galaxies. 



If we assume that the ETS are triggered by low to medium 
luminosity AGN activity and put both AGNs and ETS objects 
in a single class, one finds that the sample is largely dominated 
by galaxies with an active nucleus (78%, see Table |5]l, which 
further demonstrates that sub-mJy samples like the ATESP are 
best suited to study the evolutionary behaviour of low-power 
AGNs. 
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Fig. 12. Absolute 7?-band magnitude distribution. The sam- 
ple is divided into four different classes. From top to bottom: 
AGNs, ETS, LTS and star-burst galaxies. 



Figures 11 and 12 show the absolute magnitude distribu- 
tions in B- and R-bands for the 56 ATESP-DEEPl radio 
sources in fields DEEP la, b and c with a redshift detemiina- 
tion. AGNs are characterised by higher optical luminosities 
than ETS, LTS and SB galaxies. This is not surprising when 
we consider that in our sample all AGNs are photometrically 
and/or spectroscopically classified as quasars (see Table [6]l. 

7. Nature of the mJy and sub-mJy radio population 

In order to probe the origin (nuclear or on a larger scale) of 
the radio emission in mJy and sub-mJy sources and the phys- 
ical processes responsible for the flattening of the radio spec- 
tral index found in sub-mJy samples like the ATESP (see 
we made an overall comparison of the radio spectral index, 
the radio morphology and the optical properties of the entire 
ATESP-DEEPl sample. 

In Fig.[T3](top panel) the radio-to-optical ratio is plotted as 
a function of spectral index for the whole ATESP-DEEPl sam- 
ple (fields a, b, c an d d). The radio-to-optical ratio was defined 
foflowing lCondon | ( ,1980) , as = 5 .ioO-4('"-i2.5)^ ^here S is the 
source 1 .4 GHz flux density (in mJy) and m is the optical mag- 
nitude (here we assume the I-band magnitude). We thus can 
include sources without known redshifts. In the following we 



use both DPS and ATESP-EIS optical data (see Prandoni et al. 
2001| l, when available, while lower limits to R are given when- 
ever a source was not identified down to the limiting magnitude 
of file optical surveys (/ ~ 22.5 for EIS-WIDE and 7-24 for 
DPS DEEPl). For the sources with spectral type/redshift esti- 
mates available (either from multi-colour photometry or spec- 
troscopy) we can distinguish between ETS (red filled circles), 
LTS plus star-burst galaxies (blue stars) and AGNs (green dou- 
ble triangles). 



Figure 13 clearly shows that most of the flat-spectrum 
sources have high radio-to-optical ratios {R > 1000), typi- 
cally associated with the classical powerful radio galaxies and 
quasars. Flat-spectrum sources with low R values are preferen- 
tially identified with ETS, where the radio emission is again 
probably triggered by nuclear activity (typical radio powers 
P ~ 10^^"^^ W/Hz, see Fig. 10 and discussion therein). Star- 
forming galaxies (LTS and SB), on the other hand, are typ- 
ically associated to steep-spectrum sources, as expected for 
synchrotron emission in galactic disks or in nuclear star-bursts. 

A further radio/optical analysis of the ETS in the ATESP- 
DEEPl sample has shown that ETS with flat and/or inverted 
spectrum are preferentially compact (linear sizes d < 10-30 
kpc, see Fig. [13] middle panel). Their rather low radio lumi- 
nosities (Pi.4GHz ~ 10^^"^"^ WHz"', see Fig. 13 bottom panel) 
and the absence of emission lines in the optical spectra may 
suggest that these objects belong to the class of FRI radio 
galaxies; but FRI radio galaxies are characterised, on average, 
by steeper radio spectra and larger linear sizes (but see the lin- 
ear size - radio power relation found for B2 radio galaxies, de 
[Ruiter et al. 1 1990] and references therein). 

The compactness of the sources, together with the 
flat/inverted spectra, suggests core emission with strong syn- 
chrotron or free-free self-absorption. This could be associated 
to either very early phases of nuclear radio-activity (the 
so-called GHz peaked spectrum 



GPS - radio sources, O'Dea 



1998[ Snellen et al. 2000 1 or late phases of the evolution 
of AGNs, characterised by low accretion/radiative efficiency 
(advection-dominated accretion flow, i.e. ADAF; advection 
dominated inflow-outflow solutions, i.e. ADIOS). In the first 
case, however, larger luminosities are expected (Piaghz > 10^^ 
WHz '), while in the latter case v ery low radio po wers are 
predicted (Pjchz < 10^' Wnz '- 



see 



Doi et al. 



2005 I. Another 



intriguing possibility is that in these sources ADAF and radio 
jets coexist, as suggested for low luminosity AGNs, (LLAGNs, 
12005' and reference therein). This would 



Doi et al. 



see e.g. 

explain the somewhat brighter luminosities than expected by 
simple ADAF and can still be consistent with the presence of 



flat/inverted radio spectra (see ADAF-jet model by Falcke & 
IBierman |1999] l. 

This class of objects may also be very similar to the com- 
posite class of the so-called low power (PmuHz < 10^^^ 
WHz"') compact (< 10 kpc) - LPC - radio sources studied by 



Giroletti et al. ( 2005 1. Their host galaxies do not show signa- 
tures of strong nuclear activity in the optical (and X-ray) bands. 
Preliminary results indicate that multiple causes can produce 
LPC sources: geometrical-relativistic effects (low power BL- 
Lacertae objects), youth, instabilities in the jets, frustration by 
a denser than average ISM, and a premature end of nuclear ac- 
tivity. 

8. Summary 

In this paper we have discussed the nature of the faint, sub- 
mJy, radio population, using a sample of 1 3 1 radio sources that 
were observed at 1.4 and 5 GHz with the ATCA (the ATESP- 
DEEPl sample). A smaUer sample of 85 radio sources is cov- 



A. Mignano et al.: The ATESP 5 GHz radio survey. II. 



13 



PC 



4 - 



so 
O 




—T 1 r 1 1 r T" 



= r ^ ^ 



100 



o 



10 



T 1 1 1 1 1 1 — j — I r 



— -t--% 



_I 1 L_ 



N 

X 



26 - 



24 



- 22 



_j I |_ 



-1 

spectral index 



-f 1 1 1 1 1 1 r 



-1 r 1 1 1 [— 



Fig. 13. Radio-to-optical ratio (R, top 
panel), linear radio source size {d in 
kpc, middle panel) and 1.4 GHz ra- 
dio power (in W/Hz, bottom panel) 
against 1.4-5 GHz spectral in- 
dex for the radio sources from the 
ATESP-DEEPl sample: red filled cir- 
cles (ETS); blue stars (LTS/SB); green 
filled triangles (AGNs). Circled sym- 
bols indicate sources with extended 
and/or two-component radio morphol- 
ogy, typical of classical radio galax- 
ies. Black open circles are for iden- 
tified sources that do not have a red- 
shift/type determination. Arrows indi- 
cate upper/lower limits. Vertical dashed 
lines indicate the a - -0.5 and the 
a = values, above which source spec- 
tra are defined respectively as flat and 
inverted. Horizontal dashed lines in the 
three panels indicate, from top to bot- 
tom, values of R ^ 1000, d ^ 30 kpc 
and Pi.4GHz = 10^^ W/Hz, respectively. 



ered by deep multi-colour images. These were optically identi- 
fied down to very faint magnitudes, which was possible thanks 
to the availability of very deep multi-colour optical material (in 
U, B, V, R, I, and sometimes J and K bands). The high percent- 
age of identifications (~ 78%) makes this a sample that is well 
suited for follow up studies concerning the composition of the 
sub-mJy population and, in general, the cosmological evolu- 
tion of the various classes of objects associated with faint radio 
sources. 

We summarise our main results here. 

- For 85% of the identification sample we succeeded in de- 
riving reliable photometric redshifts, based on the available 
accurate colours (UBVRIJK). 

- Based on spectral types determined either directly from 
spectroscopy or from the photometry (or both), we find that 
at the sub-mJy level the large majority of sources are asso- 
ciated with objects that have early type (64%) and AGN 
(14%) spectra; these are of course what we would normally 
call radio galaxies and quasars. 

- Although earlier work (based on shallower optical imaging 
and spectroscopy) revealed the presence of a conspicuous 
component of late type and star-burst objects, such objects 
appear to be important only at brighter magnitudes (/ < 
19), and are rare at fainter magnitudes (19 < / < 23.5). 

- From an overall comparison of the radio spectral index with 
other radio and optical properties of the entire ATESP- 
DEEPl sample, we find that most sources with flat radio 



spectra have high radio-to-optical ratios, as expected for 
classical radio galaxies and quasars. Flat-spectrum sources 
with low radio-to-optical ratios are preferentially associ- 
ated with ETS, in which the radio emission is most plausi- 
bly triggered by nuclear activity as well, while star-forming 
galaxies are associated to steep-spectrum radio sources. 
- ETS with flat or inverted spectra are mostly compact, with 
linear size < 10-30 kpc, suggesting core-dominated ra- 
dio emission. Their low radio luminosities (in the range 
1022 _ io24 W/Hz at 1.4 GHz) and the absence of emission 
lines in their spectra (when available) suggest that they are 
FRl sources, although these would normally have steeper 
spectra and be more extended. They may therefore repre- 
sent specific phases in the life of a radio source, or may be 
similar to the low power compact radio sources discussed 
by |Giroletti et al. |p005] l. 
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Table 6. List of the optical identifications of ATESP radio sources in fields DEEPla, b and c. We list in Col. 1 the source name 
and in Cols. 2-8 the U, B, V, R, I, J, Ks-hand magnitudes (Vega system), when available. In Cols. 9 to 1 1 we list the photometric 
redshift, together with the coiTesponding SED and age (for SED templates, types and ages we refer to Hyperz documentation; 
we note that CWW and quasar templates do not provide ages). Spectroscopy information, when available, is reported in Cols. 12 
(redshift) and 13 (spectral types). References and notes are listed in Col. 14. 



Source Name 


U 


B 


V 


R 


I 


J 




Zphol 


SED 


Age (Gyr) 


Zip 


Sp. Type 


Notes 


Class 


ATESP5 


J224750-400148 


> 25.1 


21.8 


20.7 


19.3 


18.4 


- 


- 


0.35 


Burst 


8.5 


0.442 


ETS 


(a) 


ETS 


ATESP5 


J224753-400455 


> 25.1 


23.4 


22.5 


21.3 


20.1 


- 


- 


0.61 


Burst 


1.7 


- 


- 




ETS 


ATESP 


J224759-400825 


> 25.1 


25.4 


25.2 


25.2 


24.3 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J22480 1-400542 


> 25.1 


23.1 


22.1 


20.9 


19.8 


- 


- 


0.56 


Burst 


2.0 


- 


- 




ETS 


ATESP 


J224803-400513 


17.7 


18.4 


17.9 


17.6 


17.2 


- 


- 


1.00 


QSO 


- 


1.72 


QSO 


(be) 


AGN 


ATESP5 


J224809-40221 1 


22.5 


23.0 


22.7 


21.5 


20.8 


- 


- 


0.55 


Sc 


4.5 


- 


- 




LTS 


ATESP 


J22481 1-402455 


22.0 


22.7 


22.4 


21.6 


21.1 


- 


- 


0.55 


Burst 


0.05 


- 


- 




LTS 


ATESP 


J224817-400819 


>25.1 


23.9 


22.9 


21.8 


20.3 


- 


- 


0.80 


Burst 


1.0 


- 


- 




ETS 


ATESP5 


J224822-401808 


>25.I 


23.5 


22.2 


20.6 


19.4 


- 


- 


0.37 


Ell 


6.5 


- 


- 




ETS 


ATESP 


J224843-400456 


24.3 


25.6 


24.8 


24.9 


24.3 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J224850-400027 


21.8 


22.3 


21.8 


20.8 


19.7 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J224858-402708 


22.1 


22.3 


21.6 


20.5 


19.8 


- 


- 


0.48 


Burst 


0.09 


- 


- 




LTS 


ATESP 


1224911-400859 


17.8 


17.4 


16.9 


16.2 


15.6 


- 


- 


0.11 


Sbc 


- 


0.065 


LTS 


(b) 


LTS 


ATESP5 


J224919-400037 


>25.I 


22.6 


21.5 


20.0 


19.0 


- 


- 


0.35 


Burst 


10.5 


- 


- 




ETS 


ATESP5 


J224932-395801 


23.6 


24.1 


23.0 


21.6 


20.4 


- 


- 


0.60 


Burst 


2.0 


0.713 


ETS 


(ad) 


ETS 


ATESP5 


J224935-400816 


18.3 


17.6 


16.6 


15.9 


15.2 


- 


- 


O.Il 


Burst 


12.5 


0.153 


ETS 


(a) 


ETS 


ATESP5 


J224948-395918 


>25.I 


25.3 


25.9 


24.5 


23.5 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J22495 1-402035 


14.6 


18.6 


18.4 


18.0 


17.6 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J224958-395855 


13.6 


16.7 


16.4 


16.3 


15.2 


- 


- 


- 


- 


- 


0.249 


ETS 


(b) 


ETS 


ATESP5 


J225004-402412 


>25.1 


23.4 


22.3 


21.2 


20.0 


- 


- 


0.61 


Burst 


1.7 


- 


- 




ETS 


ATESP5 


J225008-400425 


18.5 


18.0 


17.1 


16.4 


15.7 


- 


- 


O.IO 


Burst 


3.5 


0.126 


ETS 


(e) 


ETS 


ATESP 


J225009-400605 


>25.I 


25.1 


23.8 


22.4 


21.0 


- 


- 


0.74 


Burst 


1.0 


- 


- 




ETS 


ATESP5 


J225028-400333 


20.7 


21.0 


20.4 


19.2 


18.5 


- 


- 


0.51 


Burst 


0.18 


0.540 


LTS 


(b) 


LTS 


ATESP5 


J225048-400147 


>24.6 


>25.7 


26.4 


25.9 


23.8 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J225056-400033 


23.9 


24.7 


24.6 


23.2 


22.2 


- 


- 


1.43 


REDQ 


- 


- 


- 




AGN 


ATESP5 


J225056-402254 


21.2 


20.4 


19.1 


18.3 


17.4 


- 


- 


0.21 


Burst 


11.5 


- 


- 




ETS 


ATESP5 


J225057-401522 


15.8 


15.3 


14.6 


14.0 


13.2 


- 


13.0 


0.01 


Burst 


0.72 


0.033 


ETS 


(af) 


ETS 


ATESP5 


J225058-401645 


>24.6 


26.3 


24.7 


23.1 


21.0 


- 


18.6 


0.96 


Burst 


1.02 


- 


- 




ETS 


ATESP5 


J225 100-400934 


>24.6 


>25.7 


26.9 


24.5 


22.5 


- 


18.0 


1.21 


Ell 


5.5 


- 


- 




ETS 


ATESP5 


J225 112-402230 


26.0 


>25.7 


27.8 


26.5 


24.4 


- 


18.5 


- 


- 


- 


- 


- 




- 


ATESP5 


J225 122-4025 24 


23.1 


23.7 


23.5 


23.1 


22.6 


- 


19.5 


2.23 


SB2 


- 


- 


- 




SB 


ATESP5 


J225138-401747 


19.4 


19.3 


18.4 


17.8 


17.0 


16.0 


14.5 


0.21 


Burst 


0.26 


0.235 


LTS 


(e) 


LTS 


ATESP 


J225206-401947 


20.3 


20.6 


20.1 


19.7 


19.1 


18.4 


17.1 


2.06 


QSO 


- 


- 


- 




AGN 


ATESP5 


J225217-402135 


22.7 


23.8 


23.4 


22.9 


22.3 


20.9 


19.2 


0.93 


BALQ 


- 


- 


- 




AGN 


ATESP5 


J225223-401841 


15.7 


15.4 


14.9 


14.3 


13.6 


12.7 


11.5 


0.04 


Sa 


12.5 


0.033 


ETS 


(af) 


ETS 


ATESP5 


J225239-401949 


14.9 


14.8 


14.3 


13.8 


13.1 


12.2 


11.3 


0.06 


Sa 


8.5 


0.033 


ETS 


(e) 


ETS 


ATESP5 


J225242-395949 


23.9 


24.2 


22.9 


21.7 


20.7 


- 


- 


0.41 


Burst 


0.36 


- 


- 




ETS 


ATESP5 


J225249-401256 


24.0 


25.7 


24.3 


22.6 


20.9 


19.4 


17.5 


0.59 


Ell 


5.5 


- 


- 




ETS 


ATESP5 


J225 3 16-40 1200 


23.7 


23.6 


22.5 


21.2 


20.0 


- 


- 


0.36 


Ell 


5.5 


- 


- 




ETS 


ATESP5 


J225321-402317 


22.1 


23.3 


22.8 


21.9 


21.2 


- 


- 


0.70 


Burst 


0.18 


- 


- 




LTS 


ATESP5 


J225322-401931 


>24.6 


25.3 


23.1 


21.7 


20.0 


- 


- 


0.30 


Burst 


10.5 


- 


- 




ETS 


ATESP5 


J225323-400453 


>24.6 


25.0 


24.0 


22.8 


21.5 


- 


18.2 


0.36 


Sa 


6.5 


- 


- 




ETS 


ATESP5 


J225325-400221 


23.1 


22.8 


21.7 


20.4 


19.1 


- 


16.1 


0.37 


Ell 


5.5 


- 


- 




ETS 


ATESP5 


J225332-402721 


>24.6 


25.6 


26.9 


24.2 


24.0 


- 


- 


- 


- 


- 


- 


- 




- 


ATESP5 


J225344-401928 


23.8 


24.9 


24.4 


24.1 


22.7 




18.8 


1.40 


Ell 


3.5 








ETS 


ATESP5 


J225345-401845 


19.9 


19.2 


18.1 


17.3 


16.3 




14.3 


0.29 


Burst 


1.7 








ETS 


ATESP 


J225351-400441 


>25.3 


23.9 


23.7 


23.2 


22.4 






2.20 


QSO 










AGN 


ATESP5 


J225353-400154 


21.3 


21.3 


21.4 


20.9 


19.7 




17.1 


2.00 


REDQ 


- 




- 




AGN 


ATESP5 


J225400-402204 


16.8 


16.4 


15.8 


15.1 


14.2 


- 


- 


0.03 


Sa 


13.5 


0.033 


LTS 


(g) 


LTS 


ATESP5 


J225404-402226 


21.4 


21.4 


20.7 


19.5 


18.6 


- 


16.2 


0.43 


Ell 


4.5 


- 


- 




ETS 


ATESP5 


J225414-400853 


15.1 


15.3 


15.0 


14.6 


14.1 


13.6 


12.8 


0.07 


Burst 


0.09 


0.032 


LTS 


(e) 


LTS 


ATESP5 


J225430-400334 


>25.3 


23.4 


22.3 


20.9 


19.4 


- 


16.4 


0.56 


Burst 


2.6 


- 


- 




ETS 


ATESP 


J225430-402329 


>25.3 


25.6 


25.1 


24.1 


22.3 




18.7 


1.03 


Ell 


4.5 








ETS 


ATESP5 


J225434-401343 


>25.3 


25.8 


25.8 


24.6 


23.2 


20.7 


19.0 


1.64 


Burst 


0.36 


- 


- 




ETS 


ATESP5 


J225436-400531 


>25.3 


22.2 


21.1 


19.7 


18.4 


17.7 


16.0 


0.44 


Burst 


4.5 








ETS 


ATESP5 


J225442-400353 


>25.3 


25.1 


24.1 


22.9 


21.0 


19.5 


17.6 


0.81 


Ell 


5.5 








ETS 


ATESP5 


J225449-400918 


23.2 


23.6 


23.3 


23.2 


22.3 


20.6 


18.9 


1.58 


Ell 


2.6 








ETS 


ATESP5 


J225509-402658 


21.6 


22.3 


21.8 


21.0 


19.4 




16.5 


0.97 


Burst 


0.13 








LTS 


ATESP5 


J2255 15-401 835 


23.3 


23.7 


23.0 


22.5 


21.1 


19.2 


15.9 


2.20 


REDQ 










AGN 


ATESP5 


J225529-401101 


>25.3 


23.3 


22.2 


21.0 


19.5 


18.7 


16.6 


0.73 


Burst 


0.36 








ETS 


ATESP5 


J224827-402515 


22.2 


22.8 


22.8 


22.0 


21.0 






1.25 


Burst 


0.36 








ETS 


ATESP5 


J225034-401936 


>24.6 


24.5 


24.0 


22.9 


21.1 






1.17 


Burst 


0.72 








ETS 


ATESP5 


J225426-402442 


>25.3 


24.1 


23.9 


23.4 


22.2 




18.7 


1.92 


REDQ 










AGN 


ATESP5 


J225443-401 147 


>25.3 


>25.9 


>25.8 


>25.7 


>23.8 


>22.2 


20.0 
















ATESP5 


J225505-401301 


20.6 


19.8 


18.6 


17.7 


16.8 


16.0 


14.7 


0.35 


Burst 


1.0 








ETS 


ATESP5 


J225511-401513 


>25.3 


>25.9 


>25.8 


>25.7 


>23.8 


21.2 


18.3 

















" On-going spectroscopy program (Prandoni et al. in prep.) 
Spectroscopy from Prandoni et al. 1 2001 1. 

Spectroscopic redshitt re-measured: published value was over-estimated. 

Same redshift obtained by Olsen et al. i 200£^: galaxy in a cluster aXz = 0.710. 

^ Spectroscopy from Vettolani et al. 1 199?^ i.^ 

^ Same redshift obtained by Di Nella et al. l| l996| . 
^ Spectroscopy obtained from fadh public data! 
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